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I. INTRODUCTION 



The flavour singlet axial current 

# = (1-1) 
is not conserved in QCD because of the triangle anomaly [1] 

d^(x)=2N f Q(x). (1.2) 
In eq. (1.2) Q(x) is the topological charge density, defined as 

Q(x) = -^^F^{x)F%{x). (1.3) 

As a consequence the corresponding Ua(X) is not a symmetry [1]. 
The non singlet partners 

JL = fyVA«V (i.4) 

are conserved, and the corresponding symmetry is spontaneously broken, the pseudoscalar 
octet being the Goldstone particles. 

If Ua(X) were a symmetry, either parity doublets should exist, or, in case of spontaneous 
breaking, the inequality < y/3m n should hold [2] . Neither of these predictions is true in 
nature, and this has been known as the Ua(X) problem for many years, before the advent of 
QCD. 

However C/a(1) is a symmetry at leading order in the expansion in ^- [3], (N c is the 
number of colours). There are arguments that the leading approximation in that expansion 
describes the main physics of QCD [4,5]. The anomaly is non leading, and can be viewed 
as a perturbation. One of its effects is to displace m rj i from zero, which corresponds to the 
Goldstone particle in the leading order approximation, by an amount which is related to the 
topological susceptibility x °f the vacuum at the leading order. The prediction is [6,7] 

2N f _ 



—X = m v + m v , - 2m K . 

J 7T 



;i.5) 



The topological susceptibility x is denned as 



X = 



J d A x(0\T(Q(x)Q(0))\0) 



(1.6) 



Leading order implies absence of fermions and in the language of the lattice this is known as 
quenched approximation. Lattice is the ideal tool to compute x from first principles, x has 
in fact been determined [8] and is consistent with the prediction of ref . [6] . An additional hint 
in favour of it is the indication that the r/' mass is higher in sectors with higher topological 
charge [9]. 

A question then arises naturally whether the Ua(X) symmetry is restored in quenched 
QCD at the same temperature at which SUa(3) is restored, i.e. at T c ~ 260 MeV [10]. 
Many models [11] predict the behaviour of the U A (l) chiral symmetry at T c . A quite general 
expectation is that the topological susceptibility should drop at T c [12], since Debye screening 
inhibits tunneling between states of different chirality and damps the density of instantons. 

Attempts have been made a few years ago to study the behaviour of x through T c [13,14]. 
The status is discussed in ref. [14]. The difficulties go back to the definition of a topological 
charge on the lattice. The correct way to define it, according to the commonly accepted 
prescriptions of field theory, is to introduce on the lattice a local operator Ql(x) for the 
topological charge density which tends to the continuum operator as the lattice spacing 
goes to zero. Ql provides a regularized version of Q(x). In going to the continuum limit 
a proper renormalization must be performed, like in any other regularization scheme. A 
specific feature of Ql is that on the lattice it is not the divergence of a current, like in the 
continuum, and hence it renormalizes multiplicatively: this means that the lattice topological 
charge of a configuration can be non integer [15]. In formulae 



As usual, P = 6/<7q- The topological susceptibility can be defined on the lattice as 



Q L = Z(f3)Qa A + 0(a 6 ). 



(1.7) 



Xl=(EQl(x)Q l (Q)). 



(1.8) 



X 
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The standard rules of renormalization then give 

XL = ^(/3)V X + M(/3) + 0(a 6 ), (1.9) 

where M(f3) is an additive renormalization containing mixings of Xl to other operators with 
the same quantum numbers and lower or equal dimensions [16]. In formulae 

M(P) = B(P)a*G 2 + P(P). (1.10) 

The terms proportional to P(@) and B(f3) are respectively the mixings to the identity 
operator and to the density of action Gi = {^iF^ v F pv ) . The additive renormalization 
comes from the singularities of the product Q(x)Q(0) as x — > and must be removed to be 
consistent with the prescription used to derive eq. (1.5) [17]. The definition of Q L is not 
unique: infinitely many operators can be defined which obey eq. (1.7) but differ by terms 
of order 0(a 6 ). The simplest definition of Ql is [18] 

Ql(x) = -^- 2 £ i^Tr (U, u (x)U pa (x)) . (1.11) 

Here e pvpa is the standard Levi-Civita tensor for positive directions while for negative ones 
the relation e^ upa = —e- pvpa holds. U pv is the plaquette in the ji — v plane. With this 
definition Z ~ 0.18 and the mixing M is large compared to the signal in the scaling region 
[19]. Z and M can be computed non-perturbatively [20-22]. Although the field-theoretic 
method is correct in principle, it is unpleasant that most of the signal is due to lattice arti- 
facts, which have then to be removed. Moreover at the time of ref. [14] the non-perturbative 
determination of Z and M was not known. 

An alternative method to determine x is the so called cooling technique [13]: the idea 
is to freeze quantum fluctuations by a local algorithm which cools the links one after the 
other. The modes relevant at a distance d are frozen after a number of steps n, which is 
proportional to d 2 , like in a diffusion process. Most of the instantons are expected to have a 
size of the order of the correlation length. After a few cooling steps, the elimination of local 
fluctuations will suppress the mixing M and make Z ~ 1, but the number of instantons will 
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be preserved, so that xl — X fl4 - m fact a plateau is reached in Ql after a few cooling steps, 
where Q L is an integer, which stands many further steps [13]. At T = and below T c the 
method works very well and agrees with the field theoretical method [16]. Approaching T c 
from below, however, the plateau becomes shorter and shorter, collapsing to a maximum, 
which at higher T decreases to non integer values. No unambiguous criterion can then be 
given to determine the value of x- The origin of the above behaviour is well understood. A 
finite temperature T on a lattice is obtained by taking a size x N T (N s spatial size, N T 
time size), with N T <ti N s . The temperature is given by 

T = 1 —. (1.12) 

As T rises to T c the correlation length becomes equal to N T a, and the instantons, being of 
the size of the correlation length, become infrared unstable. 

A third method to determine x is the so called geometrical method [23,24] by which a 
lattice configuration is interpolated by a continuous configuration on which the topological 
charge is read. In this way the charge is always an integer, which hopefully should coincide 
with the true topological charge in the continuum limit when (5 — > oo and the correlation 
length goes large. This hope, however, is not realized: at (3 — > oo lattice artefacts dominate 
[25], at least for the usual actions. No determination of x at finite T exists in the literature 
by this method. On the lattice actions can be constructed [26] which belong to the same class 
of universality of Wilson action, for which dislocations could be absent: a proof, however, 
is still lacking, that the topological charge determined by the geometrical method coincides 
with the continuum value as defined in textbook field theory. Moreover it is not clear how 
the susceptibility computed with the geometrical charge compares with the prescription of 
ref. [17] for the singularity at x — 0. Even if dislocations were absent a mixing with the gluon 
condensate, which scales with the same power up to corrections due to the dependence of 
B(f3) on (3, could still be present and should somehow be removed. Historically the method 
was proposed because eq. (1.9) had been used to estimate Xi but the factor Z 2 had been 
neglected in the analysis [18], thus leading to a value of x much smaller than what required 



byeq. (1.5). 

Recently [27] an improved operator for Q L has been defined, which has the correct 
continuum limit (1.7), but is nearer to the continuum, in the sense that Z is near to 1, and 
the mixing of xl to the identity is negligible in the scaling region. The operator has been 
tested successfully in SU(2) gauge theory with the normal Wilson action. 

What we do in the present paper is to implement the field theoretical method using 
the improved operator for Q L . Also for SU(3) we find that lattice artifacts are drastically 
reduced with respect to the original choice, eq. (1.11), and become unimportant. They are 
anyhow removed non-perturbatively. 

We redetermine x at T = 0, and study its behaviour through T c . The new determination 
of x at T = is consistent with previous determinations [19,22,29]. At finite T our main 
result is that x drops at T c by more than one order of magnitude. In sect. II the results are 
presented and discussed. The details of the determinations of Z(f3) and M(f3) are described 
in sect. III. Sect. IV contains a few concluding remarks. 

II. THE METHOD 

We will compute the lattice topological susceptibility xl defined in equation (1.9) by use 
of the improved operators defined in ref. [27] for the topological charge density 

Q?(*) = ^ E WTr (n»(x)n»(x)) . (2.1) 

In this definition, IlW is the plaquette constructed with i times smeared links Uffi(x). They 
are defined as 

U<?\x) = U lt (x), 
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The improving is a local smearing inspired by the usual cooling procedure [13]. The 
parameter c can be tuned to optimize the improvement. The formal continuum limit of 
is Q$ — °a 4 Q + 0(a 6 ) for any i. In our simulations we shall make use of the 0,1 and 2- 
improved topological charge density operators Q^\x), Ql\x) and Q^\x). Correspondingly 
we will compute the topological susceptibility 

X 

= Z^((3) 2 a\ + M«(/3) + £>(a 6 ), (2.3) 

where 

M«(/3) = B^(f3)a 4 G 2 + (2.4) 

In eq. (2.3), \ an d G 2 do not depend on the operator used for Q L ; Z^(P), B^ l \(5) 
and P^(f3) do. We checked that does not depend within errors on the lattice size, as 
expected from the fact that its value is determined predominantly by short range fluctuations 
at the UV cutoff. The multiplicative and additive renormalizations will be determined by a 
non-perturbative technique [20-22] as follows. To determine we perform a few updating 
steps on a classical configuration consisting of one instanton, and we measure Q$ at each 
step: as long as no new instantons or antiinstantons are produced and the initial instanton 
is present, the topological charge of the configuration is preserved and is equal to the known 
charge Q of the instanton. A plateau will be reached when the local fluctuations which 
produce renormalization will thermalize, and on the plateau Q$ = Z^Q, whence Z^> can 
be determined. 

To determine JtfW(/3) again we start from a configuration of known Q, e.g. the flat 
configuration (with all links U^x) = 1) where Q — 0. We then start producing from it a 
sample of configurations by the usual updating procedure used to thermalize. We measure 
at each step Xl an d as long as no instanton or antiinstanton is created or destroyed, after 
a few steps local fluctuations will be thermalized and a plateau will be reached, according 
to eq. (2.3): the first term Z^((3fa((3) A Q 2 is known (it is zero if the initial configuration is 
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flat) and can be extracted. A careful check shows that MW(/3) determined in this way 
is independent of the Q of the initial configuration: Q = and Q — 1 configurations give 
the same value of AfW(/3), within errors. In the continuum formulation of the theory this 
corresponds to the well-known fact that short distance effects, like renormalizations, do not 
depend appreciably on the classical background. 

To be sure that the number of instantons is not changed by the heating procedure, each 
configuration of the sample is checked during heating by performing a few cooling steps to 
detect its topological charge on a copy of it [28]. This is especially required at small values 
of (5 where fluctuations at the scale of the correlation length can be easily created as the 
lattice spacing is larger in physical units. In figure 1 we show the resulting topological charge 
distribution on a set of 500 initial flat configurations updated during 36 heat-bath sweeps 
at (5 — 5.90 and then frozen by 6 cooling steps. In the figure we clearly see a major peak at 
Q — and minor peaks around non-zero integer values of Q, corresponding to configurations 
in which instantons or antiinstantons have been created. The configurations representing 
these minor peaks must be discarded. To do that we operate a cut \Q\ < 5, at a S of the 
order of the width of the major peak. The background in the figure gives an estimate of the 
systematic error. We take as this systematic error the ratio A b /AQ =0 where A b is the area 
of the background in the cut interval, and Aq =0 is the area of the major peak. This error is 
added to the statistical one by quadrature. 

A similar procedure was developed for the determination of the multiplicative renormal- 
ization Z(f3), as shown in fig. 2. 

The simulations were performed on an APE QUADRICS machine and the Monte Carlo 
technique used was standard. 

For a(P) we will use the formula 



In this equation is an effective scale which in principle depends on (3 and becomes prac- 
tically constant at large enough /3, when the two-loop approximation is good for asymptotic 




(2.5) 
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scaling. As explained in next section, we will fix A^ in terms of T c and then remain in a 
small interval of (3 such that its dependence on (3 can be neglected and eq. (2.4) is valid. 

III. THE MONTE CARLO RESULTS 

We measured the topological susceptibility at zero temperature on a symmetric lattice 
16 4 and at finite temperature on a lattice 32 3 x 8. 

In figure 3 we plot the value of (xY 1 ^^ at zero temperature versus (3 for the three 
definitions Xl = 0, 1,2) of the lattice topological susceptibility. The corresponding data 
are listed in Table I. There is an excellent agreement between the three determinations. 
To fix the value of x 1 ^ 4 i n physical units we can either use [10] (3 C (N T = 8) = 6.0609(9) 
and eq. (2.4) and (1.11) which gives A L in terms of T c , or the determination of reference 
[30] A L = 4.56(11). The horizontal line is the linear fit to the 2-smear data. It gives 
(-^(V 4 ) — 175(5) MeV, and is consistent, within errors, with that of ref. [19,22,29]. The 
error includes the uncertainty in A L . 

As a check of thermalization we show in figure 4 the distribution of total topological 
charge for 5000 configurations at (3 — 6.1 for the 2-smeared charge. The figure was obtained 
by applying 6 cooling steps on fully thermalized configurations. The average topological 
charge is zero within errors, as it should be. 

In figure 5 the topological susceptibility x = (Xl ~ MW)/(ZW 2 a 4 ) at the transition 
point is shown for the 1 and 2-smeared operator, x drops by one order of magnitude from 
the confined to the deconfined phase. The results obtained with the two operators are 
compatible as they should. The data for the 0-smeared operator have very large error bars 
and are not shown in the figure. The data have been plotted versus T/T c where T c is the 
deconfining temperature. To determine T/T c we only need the ratio a((3 c )/a((3) and for that 
the two-loop expression is certainly a good approximation within the small interval of (3 
used, where A^ can be considered as a constant. 

The solid line in figure 5 corresponds to the value of x & t zero-temperature and is con- 
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sistent with the data below T c indicating that x is practically T-independent in the confined 
phase. 

In Table II we display the unrenormalized value of the topological susceptibility x^l an d 
the values for both and as a function of (3 for the measurements performed using 
the 1-smeared topological charge density . The last column is the physical value for the 
susceptibility x/^-l obtained from equation (2.3) . In Table III we show the analogous set 
of data for the 2-smeared operator. 

The quality of our 2-smeared operator can be appreciated by looking at the numbers in 
the last column of Table IV; , which includes both the mixing to the action density and 
to the identity operator (see eq. (2.4)), is 0.2-^0.3 times the subtracted signal Xl^~ . The 
deviations from constant of that ratio allows to estimate P( 2 \(3). At (3 — 6.1, P^((3)/xl^ ~ 
3%. Going to lower /3's this ratio decreases very rapidly; at larger /3's both the terms 
proportional to a 4 in eq. (2.3) die off exponentially and P < - 2 ' ) ((3) becomes dominant. A 
similar analysis on the data of Table V for the 1-smeared operator shows a much lower 
quality. For the non- improved operator the ratio M^/(x^ — M^) is much larger than 1, 
and M^{(3) is dominated by P^{(3). 

A preliminary estimate of the mixing to the gluon condensate G 2 across T c shows that 
G2 is much smoother than x &t T c and compatible with a constant. Work is in progress on 
this point to reduce the error bars. 

IV. CONCLUDING REMARKS 

We have used an improved operator for the topological charge density to determine the 
topological susceptibility of SU (3) pure gauge theory at zero temperature and its behaviour 
through T c . 

Lattice artefacts are strongly suppressed with respect to the ordinary definition, eq. 
(1.11), and are anyhow removed by non-perturbative methods. Moreover statistical fluctu- 
ations are drastically reduced. 
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Our main results are: 

1) At T = the Witten-Veneziano solution [6,7] of the Ua(X) problem is confirmed. The 
value of x is more precise than in previous determinations and agrees with them within 
errors. 

2) x drops to zero at T c . 

The method used rests on basic concepts in field theory, and on the assumption that 
lattice is a legitimate regularization of it. No additional assumptions are needed. 
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Figure captions 

Figure 1. Distribution of topological charge Q for a set of 500 configurations obtained by 36 
heat-bath updatings of the flat configuration and 6 cooling sweeps. f3 = 5.90. 

Figure 2. Distribution of topological charge Q for a set of 2000 configurations obtained by 15 
heat-bath updatings of a 1-instanton configuration and 6 cooling sweeps. f3 = 5.75. 

Figure 3. x & t T = 0. The straight line is the result of the linear fit of the 2-smeared data. The 
improvement from to is clearly visible. 

Figure 4. Topological charge distribution for an ensemble of 5000 fully thermalized configurations 
at P = 6.1, (2-smeared operator). 

Figure 5. versus T/T c across the deconfining phase transition. The horizontal band is the 

determination at T = of figure 2. 

Table captions 

Table I. x 1 / 4 y 1 Al from the 0,1 and 2-smeared operators on a 16 4 lattice. 

Table II. T/T c , X l\ M« and as a function of (3 for the 1-smeared operator. 
Table III. T/T c , Xl\ and as a function of (3 for the 2-smeared operator. 
Table IV. Data for Xl? ano ^ for the 2-smeared operator. 

Table V. Data for Xl 1 an d for the 1-smeared operator. 
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Table I 



p 


V (V4) /* 
AO— smear/ ii i 


v (l/4) /a 
Al— smear/ ii i 


v (l/4) /a 
A 2— smear / -"VL 


5.90 


36.7(12.9) 


39.1(1.1) 


38.8(0.8) 


6.00 


40.1(14.8) 


39.1(1.1) 


38.4(0.8) 


6.10 


43.1(14.4) 


39.0(1.2) 


38.1(0.8) 



Table II 



p 


T/T c 


10 5 x 


10 5 x ilfW 


Z 0) 


10 6 x X /A| 


5.90 


0.834 


2.35(7) 


0.88(6) 


0.36(2) 


2.46(32) 


6.00 


0.934 


1.55(5) 


0.62(3) 


0.39(2) 


2.07(25) 


6.06 


0.999 


1.05(3) 


0.53(4) 


0.40(2) 


1.38(19) 


6.10 


1.045 


0.728(21) 


0.48(3) 


0.42(2) 


0.76(13) 


6.22 


1.197 


0.428(12) 


0.34(2) 


0.45(2) 


0.38(11) 


6.30 


1.310 


0.334(9) 


0.314(15) 


0.48(2) 


0.11(10) 


6.36 


1.402 


0.293(10) 


0.281(10) 


0.49(2) 


0.08(10) 
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Table III 



p 


T/T c 


10 5 x x f 


10 5 x M( 2 ) 




10 6 x X /A| 


5.90 


0.834 


3.27(9) 


0.71(6) 


0.49(2) 


2.07(18) 


6.00 


0.934 


2.04(7) 


0.46(3) 


0.51(2) 


2.00(18) 


6.06 


0.999 


1.31(4) 


0.38(3) 


0.53(2) 


1.54(15) 


6.10 


1.045 


0.76(2) 


0.33(2) 


0.55(2) 


0.85(9) 


6.22 


1.197 


0.320(9) 


0.227(14) 


0.58(2) 


0.32(6) 


6.30 


1.310 


0.227(6) 


0.197(9) 


0.60(2) 


0.10(4) 


6.36 


1.402 


0.184(6) 


0.168(7) 


0.622(13) 


0.07(4) 



Table IV 



p 


10 5 x X ( L ] 


10 5 x 




5.90 


3.51(7) 


0.71(6) 


0.253(23) 


6.00 


2.18(5) 


0.46(3) 


0.270(20) 


6.10 


1.39(3) 


0.33(2) 


0.315(22) 



Table V 



p 


10 5 x X ( L ] 


10 5 x M« 




5.90 


2.48(5) 


0.87(4) 


0.540(32) 


6.00 


1.64(4) 


0.60(2) 


0.576(29) 


6.10 


1.12(2) 


0.474(13) 


0.734(34) 



16 



-a 




on; 
c 

i-S 

CD 



I — I 

o 

H 



200 - 



> 

« 180 V 



160 - 



140 ■ 



♦ smear 

• 1 smear 
■ 2 smear 




- 45.5 



CD 



- 36.4" 



- 31.8 



5.85 



5.90 



5.95 



6.00 6.05 



P 



6.70 6.75 



APPENDIX 
Erratum to: 

"Topological susceptibility at zero and finite T 
in SU(3) Yang-Mills theory" 

[Nuclear Physics B494 (1997) pg. 281] 

B. Alles, M. D'Elia, A. Di Giacomo 

In the computation of Table 1 and of the last columns of Tables 2 and 3 some trivial 
numerical error was made. These numbers are computed from the lattice data of Tables 2-5 
by use of Eq.(l.lO) and (2.5). The correct tables are displayed below showing a marginal 
violation of asymptotic scaling. 

It is therefore better to present our results in terms of the scaling ratio x^V^c, discarding 
the approximation Eq. (2.5) and using instead the numerical determinations of T c at different 
(3 of Ref. [1]. The results are shown in Figure 3 and 5 below. The conclusions of the paper 
stay unaltered. 

The MeV units in Fig. 3 have been obtained by using the ratio T c /y/a = 0.629(3) 
from Ref. [1] and the value y/a = 420 MeV. From the 2-smeared data we derive x 1 ^ 4 = 
170(7) MeV. If the ratio T c /y/a is taken from [2], this result becomes x 1/4 = 174(7) MeV. 
The new Figure 5 displays the scaling ratio xl^t- The temperature axis in this Figure has 
been determined again by using T c from Ref. [1]. 



Table 1 



p 


v (l/4) /* 
AO— smear/ ii £ 


v (l/4) /* 
Al— smear/ ii £ 


V (V4) /* 
A 2— smear / ii i 


5.90 


42.3(10.4) 


40.2(1.2) 


39.6(0.9) 


6.00 


37.8(8.2) 


38.8(1.1) 


38.4(0.8) 


6.10 


40.6(7.8) 


37.1(1.0) 


36.7(0.7) 
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First and last columns in Table 2 



o 

13 


10 x x/Al 


5.90 


2.39(30) 


6.00 


2.02(24) 


6.06 


1.41(20) 


6.10 


0.73(13) 


6.22 


0.39(11) 


6.30 


0.11(10) 


6.36 


0.08(10) 



First and last columns in Table 3 



p 


10- 6 x X /Ai 


5.90 


2.24(21) 


6.00 


2.00(18) 


6.06 


1.43(13) 


6.10 


0.74(7) 


6.22 


0.25(5) 


6.30 


0.11(4) 


6.36 


0.07(4) 
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Figure 3 



211.3 



^ 184.9 



X 



158.5 



132.1 



♦ O-smear 

• 1 -smear 
■ 2-smear 



0.8 



0.7 



0.6 X 
0.5 



5.85 5.90 5.95 6.00 6.05 6.10 6.15 



P 



Figure 5 



0.20 



0.15 



0.10 



0.05 



- x/t; t=o 

• l-smear 
■ 2-smear 



1 



0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

T/T c 

We are grateful to Ettore Vicari for pointing out to us the inexact numbers in Table 1. 
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